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Quantum Phase Transition from Superfluid to Localized Condensates of Bose Fluid in 
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We develop a model of a strongly correlated Bose fluid model in a confined potential for the 
purpose of analyzing the localization of Bose-Einstein condensation and the disappearance of super- 
fluidity. This work is motivated by the recent observation of a quantum phase transition in liquid 
*He in porous glass at high pressures. By introducing a pressure-dependent localization length of 
the localized condensates, we could make a new analytical criterion for the localization of the con- 
densate. Without introducing any free parameters, the resulting critical pressure of the transition 
from superfluid to localization is shown to be quantitatively consistent with observations. 

PACS numbers: 67.40.-w, 05.30.Jp, 64.60.Cn 



CZ3 



I 

o 
o 



> 

cn 
m 
o 

o 



o 
o 



X 
J3 



Bose-Einstein condensation (BEC) and superfluidity of 
liquid ^He in random environments such as porous glass 
are active problems in quantum fluid research. In par- 
ticular, finding out how disorder competes against long- 
range order of a Bose fluid has stimulated both experi- 
mental and theoretical studies. 

Below 2.17 K, Hquid ^He enters both BEC and in- 
viscid superfluid states. Superfluidity is a macroscopic 
quantum phenomenon like superconductivity, and under- 
standing them has been a major goal of quantum statisti- 
cal physics and low temperature physics. Various distinct 
features of superfluidity have been explained using the 
phenomenological two fluid model in which the system 
consists of the inviscid superfluid and a viscous normal 
fluid. On the other hand, the BEC state of liquid *He, in 
which a macroscopic number of particles occupies a single 
particle ground state and is described by a macroscopic 
wave function, was confirmed by neutron scattering ex- 
periments 12] ■ The inviscid superfiow can be understood 
using this BEC wave function However, the relation 
between the BEC and superfluidity is not necessarily un- 
derstood; the one is neither necessary nor sufficient for 
the other. For example, in a two-dimensional Bose sys- 
tem, Kosterlitz and Thouless showed that superfluidity 
can exist even without a BEC Q , and this state was later 
observed in ''He films 0. 

The random environment of a porous glass is a good 
system for studying the relation between BEC and su- 
perfiuidity ^] . In general, BEC in a random system has 
received considerable attention because localization ef- 
fects allow some condensed particles to join the normal 
fluid rather than the superfluid. The phase diagram of 
this system has been discussed, showing a specific non- 
superfluid phase 0. Porous glass is often used as a 
random media in experimental studies. Vycor glass or 
Geltech silica is porous, having 30 to 70 % of its vol- 
ume containing wormholelike pores with diameters rang- 
ing from 25 to 100 A. By adjusting the pore size, one 
can change the strength and characteristic scale of ran- 
domness, whereas adjusting the adsorbed ''He coverage 
or the applied pressure effectively changes the density of 
^He from the weakly to the strongly correlated region. 



By using torsional oscillators, Reppy et al. showed some 
remarkable features of the superfluid component and the 
superfluid critical temperature for changing pore sizes of 
porous glass or coverage 0, |^ • The superfluid density in 
such porous glass is smaller than that of bulk ''He and its 
critical temperature decreases with the coverage. Below 
a certain coverage, the superfluid density can no longer 
exist, even near K. These results show that superfluid- 
ity is broken by the random environment. 

It is also important to find out how disorder affects the 
BEC. BEC and its elementary excitation in liquid ''He 
can be observed by inelastic neutron scattering. Bulk 
^He has typical excitations consisting of phonons and 
rotons 0. Dimeo et al, Plantevin et al. and Pearce 
et al. used neutron scattering and a torsional oscillator 
to observe the elementary excitations and the superfluid 
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transition, respectively, of He in porous glass 
Surprisingly, Plantevin et al. found that these elemen- 
tary excitations existed even above the superfluid criti- 
cal temperature, which is a clear evidence of BEC with- 
out superfluidity. Recently, Yamamoto et al. observed 
the disappearance of superfluidity of liquid "^He confined 
in porous Geltech silica of pore size 25 A at pressures 
P > 3.5 MPa and temperatures near K without solidi- 
fication 12], suggesting a new quantum phase transition 
to localized BECs under the effect of a strong correlation. 
For Vycor glass with larger pores, this disappearance of 
superfluidity is suppressed by the solidiflcation of liquid 
^He 

A Bose system in a random environment is also an in- 
teresting system to study theoretically. The long-range 
order correlation due to BEC can compete with the dis- 
order, so that the BEC and its superfluidity are sup- 
pressed. Huang and Meng proposed a model for a three- 
dimensional dilute Bose gas in a random potential con- 
sisting of a small coverage of *He in porous glass [isj . 
By introducing the size dependence of the random po- 
tential into their model, we could quantitatively explain 
the experimental results 7] observed by Reppy et al; 
these results included the disappearance of superfluidity, 
even at the finite coverage and zero temperature [ij i and 
the presence of BEC below the critical coverage of su- 
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perfluidity. Another model is the Bose-Hubbard model. 
This model is effective in the strongly correlated region 
at small densities; it includes the transfer energy, a ran- 
domly distributed on-site energy, and the on-site repul- 
sion. Fisher et al. showed that the Bosc-Hubbard model 
predicts the existence of the Bose glass phase, the su- 
perfluid phase, and the Mott insulating phase The 
Bose glass phase is similar to the Anderson insulating 
phase in metals In the Bose glass phase, the con- 

densed particles are localized and thus do not contribute 
to superfluidity. Although some theoretical studies 
examined the Bose glass phase, there is no clear experi- 
mental evidence for this phase. 

This system has also been studied using a trapped 
alkali-atomic BEC confined in optical lattices and optical 
speckle beams 18, 19] . The Bose-Hubbard model and 
the Gross-Pitaevskii model have been used to advance 
our understanding of this system [l^il^ll^; however, 
some features are incompletely understood, particularly, 
the superfluidity and the Bose-glass phase. 

In this paper, we focus on Yamamoto et aVs ^12] obser- 
vation of the disappearance of superfluidity at high pres- 
sures and very low temperatures. This disappearance of 
superfluidity can be considered to be a pure quantum 
phase transition that is induced by the strong correla- 
tions and the confinement under the conditions of the 
experiment. However the Bose-Hubbard model is not a 
good model to use for these conditions for the following 
two reasons. First, this transition is very sensitive to pore 
size and it is suppressed by solidification of liquid '*He for 
larger pores. We thus have to consider the pore size de- 
pendence, which is lacking in the Bose-Hubbard model. 
Secondly, we have to consider the liquid phase at high 
pressures - conditions under which the Bose-Hubbard 
model fails because it has only a single energy level per 
one site. 

We start from a model of three-dimensional Bose fluid 
in a confined potential that depends on pore size. Then 
we introduce a localization length Lg for localized BECs 
and propose a new analytical criterion for the system. 
We calculate the energy of the localized BECs as a func- 
tion of Lg, and find the value Lg = L™'" that minimizes 
the energy. The elevated pressure reduces i™™ from the 
system size down to the order of the pore size. Then 
we determine the critical pressure above which be- 
comes comparable to the pore size and the superfluid 
component disappears. The critical pressure is found to 
be quantitatively consistent with the experimental value 
[T^ . without introducing any free parameters. This re- 
sult gives strong evidence that the observed disappear- 
ance of superfluidity is caused by the transition to local- 
ized BECs. 



Bose-fluid in the confined potential is 

H-nN = ^[£(fe) - ^]aHk)a{k) 
k 



i U{k^-k2)a\ki)a{k2) 



+ ^ a^ki + q)aHk2 - q) 

ki ,/u2 ,Q 

X go{q)a{k2)a{ki), 



(1) 



where a{k) and at(fe) are respectively the free parti- 
cle annihilation and creation operators with the discrete 
wave number k = 2nn/L characterized by the integer 
n — {nx,ny,nz). The kinetic energy is e(fc) = fi.^/c^/2m 
with m the mass of one particle, the chemical poten- 
tial, U{k) the external confined potential, go{k) the in- 
terparticle interaction, and V = the volume of the 
cubic system with size L. We assume that number Nq 
of particles are condensed at the lowest wave number 
fc = and thus replace the operators of Bose particle 
a{k — 0) and a'^{k ~ 0) with VIVq. The second and third 
terms of the right-hand side of Eq. |^ represent the in- 
teraction between a particle and the confined potential 
and the interparticle scattering, respectively. Here we 
treat U{k) using the second-order perturbation theory 
and goik) by the ring-approximation, which is one of the 
simplest method for strongly correlated system of high 
density. For these calculations, we introduce the Green 
function 



G(fco,fe) 



Gii(fco,fc) Gi2(A:o,fe) 
G2i(fco,fc) G22{kQ,k) 



(2) 



with k 0. Calculation of G starts from the non- 
perturbative Green function 



GO(fco,fe) 



G°2iiko,k) GUko,k) 

n/ako,k) 

h/^{-ko, -k) 



(3) 



with ^(fco, fe) = hko — e{k) + ^ + \0. The effective inter- 
particle interaction g{qo, q) with the ring approximation 
is given by 



5(90, g) = 



9o{q) 



l-^g,{q)GlM.Q) 



(4) 



The Green function &(kQ,k) that includes this effective 
interparticle interaction can be obtained from the follow- 
ing Dyson equation: 



The grand canonical Hamiltonian of three-dimensional G^iko, fe) = G°(fco, k) + G°(fco, fe)S(fco, fc)G^(fco, fe), (5) 
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where S(fco, k) is the self energy defined as 
t{ko,k) 



S2l(fc0,fc) S22(fc0,fc) 

I]ii(fco,A;) =i;22(-fco, -fc) 
No 



(6a) 



[5o(0)+go(fc)] 



X G?i((7o,q)g(fco -qo,k~ g), 



Si2(fco,fc) =S2i(fco,fc) = —go{k). 



(6b) 
(6c) 



FinaUy, we can obtain the Green function G'ii(fco, fc) with 
the second-order perturbation of U{k) as foUows: 

Gn{ko,k)=Gf,\ko,k) + Gf?{ko,k), 



Gfi{ko, k) 
Gf^{ko,k) 



{27ri)S{ko)No 



[G\,{ko,k)+G\^{ko,k)] 



X [GUko,k) + Gl,{ko,k)]\Uik)\', 
Y,[Gniko,k)GUko,q)G\,{ko,k) 



1 



+ G\,iko,k)G{,{ko,q)Gl,{ko,k) 
+ G\^iko,k)Gl^iko,q)G{^{ko,k) 
+ G\^{ko,k)Gl2{ko,q)Gl,{ko,k)] 

x|{7(fc-q)p. (7) 

The chemical potential is given by the Hugenholtz-Pines 
theorem p^ : 

N Nn 



TFTr^Jl dqo G?l((?o,g).9(-'7o,-g)■ 
l/ iTTiyq^Q 



(8) 



The number of condensate particles and the total energy 
are 



No + dko Gniko, fc)e+'""°, 



(9) 



- dfco [infco + £(fe)+^]Gii(fco,fc)e+'^°°. 

(10) 

Here N is the total number of particles. 

The calculation of superfluidity is based on the two 
fluid model in which the number N of particles equals 
the number in the normal fluid component Nn and the 
superfluid component Ng. The superfluid component Ng 
can be calculated using linear response theory Be- 
cause of its viscosity, only the normal fluid responds to a 
small, applied velocity field. Thus, the normal fluid com- 
ponent can be defined by the response of the momentum 



density to the external velocity field. After some tedious 
calculations, the linear response theory gives the final 
form of Nn and Ng as 

X [Gii(fco, fc)2 - Gi2(fco, fc)G2i(fco, fe)]e+"^-°", 

Ng^N-Nn. (11) 

In the derivation, we assumed the rotational invariance. 
This derivation of Nn and Ng is consistent with the Kha- 
latnikov's method based on Galilean invariance |2||. 

Introducing a new criterion, we consider localized 
BECs. We divide the total cubic system of the volume 
V ^ into cubic subsystems of the same shape except 
the volume Vg = L^. Then we assume that every cubic 
subsystem has the same localized BEC and the number 
of particles Ng = NVg/V. The energy Eg of one subsys- 
tem can be obtained by replacing the volume V and the 
number of particles N in the above equation with Vg and 
A^g. Because the number of the localized BECs is V/Vg, 
the energy of the total system should be EgV/Vg. The 
ideal or weakly interacting Bose gas has Eg that is pro- 
portional to Vg, whereas in the strongly correlated Bose 
fluid. Eg is no longer a simple linear function of Vg-, thus, 
Eg/Vg depends on Vg. Calculating Eg/Vg as a function 
of , we determine the value Vg — V^™™ that minimizes 
the total energy E = EgV/Vg. Then, the localization of 
BEC can be defined by the following criterion. 

• If yj"™ exceeds V, the system is a non- localized 
BEG state. 

• When Vg'"'" is reduced to the pore size, we assume 
that the BEG is local ized within the localization 
length Lf'' = ^Fg""" 

To apply this criterion to the experiment by Yamamoto 
et al. [12j , we numerically calculated the condensate par- 
ticle A'o, superfiuid component Ng and the localization 
volume FJ"'" as a function of the pressure P = —dE/dV 
by changing the number of particles TV, and then com- 
pared our results with Yamamoto et aVs experiment 
|l2|. To make a quantitative comparison, we used the 
following numerical parameters. The mass of a ^He 
atom is TO ~ 6.6 x 10^^^ kg. For the bare interpar- 
ticle interaction, we used the simple gaussian interac- 
tion go{k) = voicV^n)^ exp[— fc^(T^/2], where parameters 
vo — 2.1 X 10^^'' J and cr ~ 2.5 x 10~^^ m were deter- 
mined using a comparison to the effective interparticle 
potential |2y| proposed by Aziz et al. The total volume 
of Geltech silica used in Yamamoto et aVs experiment 
was V ~ 5.9 X 10~®m'^. For the external confined poten- 
tial U{k), we used the random distribution 



{U*{k)U{k)) = C/^exp 



47r2 



(12) 



We then calculate physical values for many ensemble of 
U{k) and took an arithmetic average. In Eq. H12|l . 
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FIG. 1: Calculation results, (a) Dependence of the conden- 
sate particle number A'^o, the superfluid component A's, and 
the volume V^""° of the localized BECs on the pressure, (b) 
Log-plot of the dependence of L^'" • on the pressure 



Tp ~ 2.5 X 10"^ m is the pore size of Geltech silica. 
We chose the characteristic strength of the potential 
Uq — 2.0 X 10~^^ J to ensure the critical adsorbed liq- 
uid "^He coverage ric ~ 20/imol/m^ per pore surface area 
below which superfluidity also disappears, as we did in 
our previous study 

Because all parameters were fixed, we can quantita- 
tively compare calculations to experiment. Figure ^a) 
shows the pressure dependence of iVo, A's, and Vg™'". At 
the pressure P ^ 4.2 MPa, the superfluid component iVs 
disappears and V™™ decreases dramatically. We there- 
fore define this pressure as the critical pressure Pc- The 
condensate particle number Nq also decreases just below 
Pc- This number recovers above P^, although the con- 



densates are localized. The dependence of i™'" in Fig. 
n (b) shows that the locahzation length reduces to 
about the pore size Vp. Moreover, the critical pressure 
Pc ^ 4.2 MPa is nearly equal to the experimental value 
of Pc ~ 3.5 MPa. Therefore, we conclude that the exper- 
imental disappearance of superfluidity at high pressures 
is caused by the transition from a normal BEG to local- 
ized BECs. We made a similar calculation for the case of 
porous Vycor glass, which has larger pores than Geltech 
silica. Using the numerical parameters Tp ~ 7.0 x 10^^ m 
and Uq ~ 8.9 x 10^'^^ J, we obtained the much larger 
critical pressure of Pc ~ 9 MPa, which is too high for a 
BEG to be localized against solidification of liquid *He. 
This is also consistent with experimental results. 

In conclusion, we investigated the localization of BEG 
and the disappearance of superfluidity in a strongly corre- 
lated Bose fluid at high pressures. In the model, the BEG 
is assumed to be localized in the pores when the pressure- 
dependent localization length equals the pore size. With- 
out introducing any free parameters, we showed that the 
critical pressure above which BEG localizes is quanti- 
tatively consistent with the experimental result of liq- 
uid ^He in porous Geltech silica. Our criterion can be 
generally applicable to many systems with localization 
or quasi-condensation. For example, this formulation 
may be closely connected with the renormalization group 
analysis of the Kosterlitz-Thouless transition Q ■ 

We acknowledge Keiya Shirahama for many useful dis- 
cussions. 
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